Magnetic fields from inflation: the transition to the radiation era 
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We compute the contribution to the scalar metric perturbations from large-scale magnetic fields 
which are generated during inflation. We show that apart from the usual passive and compensated 
modes, the magnetic fields also contribute to the constant mode from inflation. This is different 
from the causal (post-inflationary) generation of magnetic fields where such a mode is absent and 
it might lead to significant. non-Gaussian CMB anisotropies. 
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I. INTRODUCTION 

Magnetic fields are observed in cosmic structures over 
a wide range of length-scales and redshifts; from galaxies 
to regions around high rcdshift quasars, from clusters and 
superclusters to low density filamentary regions [l[ . The 
field values are a few microgauss in galaxies and clusters, 
and of the order of the nanogauss in filaments. Recently, 
Fermi and HESS data have been used to put a lower 
bound of at least 10 _17 G on the intensity of magnetic 
fields in the intergalactic medium and even in voids Q. 
Finding an explanation for these magnetic fields is chal- 
lenging, and their origin to date remains an open prob- 
lem. 

One possibility is that magnetic fields have been gener- 
ated in the primordial Universe 0] . In particular, primor- 
dial magnetogenesis mechanisms operating during infla- 
tion have the advantage to provide magnetic seeds filling 
the entire Universe, possibly with significant amplitude 
also at very large scales. This goes in the right direction 
to explain both the ubiquity of the observed fields and 
the uniformity of the measured amplitudes. 

In this paper, we focus on inflationary magnetic 
fields, generated by breaking conformal invariance of 
electromagnetism via a term in the action of the form 
f 2 (<p)F fJ,,y F^ v . This coupling of the electromagnetic field 
to the inflaton was first proposed in Refs. [J, Q, and 
subsequently reanalyzed in [6[ for different categories of 
string-inspired inflationary scenarios (see also [7|)- The 
same kind of action has been considered in the context 
of dilaton electromagnetism M and DBI (Dirac-Born- 
Infcld) inflation Q. In Ref. [l(| it has been pointed out 
that, since f(<p) plays the role of the inverse coupling 
constant to the charged Dirac field (e.g. the electron), 
it must remain large because perturbation theory of the 
interaction of the Dirac field with the electromagnetic 



field is only trustable in the small coupling regime. This 
constraint greatly reduces the capability of the model to 
give rise to a significant magnetic field amplitude. Re- 
cently, it has been pointed out that this problem could 
be circumvented by coupling not only the electromag- 
netic field but the entire matter Lagrangian to the in- 
flaton [lTj |. However, a multiplication of the entire mat- 
ter Langrangian with f 2 (ip) could be absorbed in a field 
redefinition and will therefore not lead to any physical 
effects. Another possibility to circumvent the problem 
would be to multiply the coupling term with f(ip), i.e. 
y^ej^A^ip — > y 7 jef(ip)'f fJ 'A tI ip. This ensures that the cou- 
pling to the canonically normalized field, f(ip)A^ remains 
constant, however this coupling explicitly breaks gauge 
invariance which is only recovered when f(<p) freezes in 
after inflation. 

In this work we still concentrate on a coupling of the 
form / 2 (</j)F /il 'F Alzy , because it has the advantage to be 
quite general and simple, and because we believe that the 
main feature of the result we obtain does not depend on 
the specific form of the coupling. Note that helical mag- 
netic fields can be generated by coupling a pseudoscalar 
inflaton to the FF-term. However, it has been shown 
that this coupling generically leads to blue spectra, which 
do not have enough power on large scales to be the seeds 
of the large-scale coherent fields observed in galaxies and 
clusters |12| . 

Starting from the consistent assumption that the elec- 
tromagnetic field arising from the amplification of vac- 
uum fluctuations is subdominant and does not affect 
the background dynamics of inflation, previous analy- 
ses have evaluated the spectrum of the electromagnetic 
energy density and studied the conditions under which 
this kind of coupling gives rise to interesting magnetic 
field amplitudes after inflation 0,0. However, the con- 
tribution of the electromagnetic energy density, even if it 
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does not affect the background, still affects metric per- 
turbations at first order in perturbation theory. There- 
fore, here we proceed one step further and calculate the 
scalar metric perturbations induced by the electromag- 
netic field at first order in perturbation theory on su- 
pcrhorizon scales. We assume that the electromagnetic 
energy-momentum tensor is first order, and consequently 
the electromagnetic field is half order. Another possibil- 
ity would be to set the electromagnetic field first order, 
and its energy-momentum tensor second order. However, 
the metric perturbations induced by the electromagnetic 
field turn out to be much larger than second-order pertur- 
bations in the inflaton; the latter can therefore be consis- 
tently neglected, effectively going back to the first-order 
scheme. The point of view of treating the electromag- 
netic energy-momentum tensor as a second-order pertur- 
bation has been taken in Ref. [l3[ where, nevertheless, 
the second-order inflaton perturbations have also been 
neglected: consequently, Ref. [HI effectively adopts the 
perturbation expansion used in this work. 

In the present analysis we find that during inflation, 
the large-scale solution for the Bardeen potentials is 
, $_ ~ Q^(k,r])/(kr]) 2 , sourced by £2^ which de- 
notes the electromagnetic anisotropic stress normalized 
to the background energy density (note that here we iden- 
tify the metric perturbation with their r.m.s. (root of 
the mean square) amplitude, c.f. discussion in Section 
IIV A[) . Therefore, even though fl^j <C 1, the Bardeen po- 
tentials become large on supcrhorizon scales. However, 
we demonstrate that the ratio of the Weyl to Ricci ten- 
sors is small, because it is determined by fij^ -C 1. It also 
turns out that perturbations in other gauges, like comov- 
ing and synchronous gauge do remain small: therefore 
perturbation theory is valid during inflation. 

We then match the inflationary solution for the metric 
perturbations at the end of inflation to the solution in the 
radiation era on the surface of constant background en- 
ergy density, in the usual way [l4j ]. From this, we derive 
the metric perturbations at superhorizon scales during 
the radiation-dominated era. The matching shows that 
the large, l/(fc?y*) 2 contribution to the Bardeen potential 
is transferred entirely to the decaying mode in the radia- 
tion era. On the other hand, \&+ during the radiation era 
gets a contribution at next order in the large-scale expan- 
sion, i.e. at order 0((krj*) ). This constant term adds 
to the usual 'passive' and 'compensated' modes leading 
to a new effect in cosmic microwave background (CMB) 
anisotropics, and it is specific to inflationary generated 
magnetic fields (it is absent if the magnetic field is gen- 
erated by a causal process). 

The rest of the paper is organized as follows: in Sec- 
tion |TT] we present the perturbed Einstein equations and 
derive the Bardeen equation in the presence of a nonzero 
electromagnetic source; in Section lPHl wc assume a power- 
law evolution in time for f(ip), and calculate the scalar 
electromagnetic anisotropic stress which represents the 
dominant source term of the Bardeen equation at su- 



pcrhorizon scales; in Section IIV1 we find the solutions 
for the Bardeen potentials both in the inflationary and 
radiation-dominated eras, and perform the matching. In 
Section IIV Dl we derive the solutions at next-to-leading 
order in the large-scale expansion and in Section [V] we 
conclude. Some details of the calculations are deferred 
to appendices. 

Notation: Throughout this paper we use conformal 
time ?7, comoving space coordinates x and wave vectors k 
with the metric ds 2 = a 2 (r))(—drf -\-5ijdx i dx^)\ 4d space- 
time indices are Greek letters while 3d spatial indices are 
Latin letters and spatial vectors are denoted in bold face. 
For the metric and scalar field perturbations we follow 
the conventions of while for the electromagnetic ac- 
tion and field quantization we follow the conventions of 
[flUBl- We define the Planck mass by rap = (\/87rG) _1 . 



II. METRIC PERTURBATIONS SOURCED BY 
THE ELECTROMAGNETIC FIELD 

We consider an electromagnetic field generated during 
inflation by breaking of conformal invariance, as speci- 
fied in Section Hill We assume that the electromagnetic 
energy-momentum tensor is first order in perturbation 
theory, meaning that the electric and magnetic fields are 
half order (c.f. discussion in the Introduction). Dur- 
ing inflation, the background evolution is therefore de- 
termined only by the background scalar field, whereas 
both the scalar field perturbation and the electromag- 
netic field contribute to the first order energy-momentum 
tensor, ST a p = 5T v a p + T° m/3 . The scalar field driving 
inflation is split into a background part and a first-order 
perturbation as y(x, rj) = <£n(?7) + ^¥ , ( x ; r ?) i so that the 
background equations are [15[ 

47rG(^) 2 = n 2 -W and ip%+2H<p' +a 2 V v> = . (I) 

A prime denotes derivative with respect to conformal 
time r\ and 1-L = a' /a. Since the electromagnetic field is 
half order, the electromagnetic energy-momentum tensor 
can be decomposed with respect to the unperturbed ve- 
locity of the FL background, u a = a _1 (l, 0), and to the 
unperturbed metric g a p , and it is gauge invariant fl7l ] : 

= (Pc B ,+Pc ia )u a uP+p Bm g ap + 2u^l + n°£ . (2) 

The parentheses around the superscripts of the third 
term indicate symmetrization. We focus on scalar pertur- 
bations and use longitudinal gauge with the notation (l5| 

ds 2 = a 2 [-(1 + 2$)d?7 2 + (1 - 2*)dx 2 ] . (3) 

Note that the 'names' <J> and \& are interchanged with re- 
spect to The perturbed Einstein equations, SG a g = 
8TrGST a g , for scalar perturbations in Fourier space in the 
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presence of an electromagnetic field arc: 

3W + {2V 2 + H')$ + k 2 ^ = 
- AirG{ip' 8(p' + V v a 2 5(pj - 4irGa 2 p cm , 

k 2 

+ 2ft*' + ft$' + (2ft 2 + U')<$> - — ($ - *) 

o 

AirGup' Q 8<p' - V^a 2 5(p*j + AnGa 2 p cm , 

k j 

+ ft* = 4>irG<p' 5<p - 47rGai — g omj - , 
fc 2 ($ - *) = -87rGa 2 n s . 



(4) 



(5) 

(6) 
(7) 



Here p em (k), p e m(k) and q em j(k.) are the electromag- 
netic field energy density, pressure and Poynting vec- 
tor in Fourier space obtained from Eq. ©. Ils(k) = 
— 3/2k l kjH c J nli ('k) is the scalar part of the electromag- 
netic anisotropic stress; it is of the same order of mag- 
nitude as the electromagnetic energy density. These 
equations can be combined into a second-order evolution 
equation for the variable "J, the Bardecn equation, 



2ft'- 



+ k 2 )^ = S 



(8) 

The source term S is due to the presence of the electro- 
magnetic field. It is given by 



S = 8nGa 2 



a 2 k 2 



2 {Pern Pern) 
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k 2 
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<P'o) 


k 2 a 







(9) 



We want to solve Eq. © , in order to determine the effect 
of the electromagnetic field on the scalar metric pertur- 
bations. We are interested in the solution at very large 
scales k\-q\ <SL 1. We consider slow-roll inflation with 



a ~ a\ 



l+e 



n 



n 



(10) 



at first order in slow-roll. The slow-roll parameters e and 
e 2 are defined by [l8| 



H 2 - ft' = eV. 2 , 
From this we infer 



e' = 2e(3e 2 + 2e)U 



ft = (3e 2 + e)ft. 



(11) 



(12) 



Using these expressions together with Eqs. (JTJ), and defin- 
ing the new variable x = \krj\ = —krj, Eq. ([SJ can be 
rewritten as 



d 2 ^ 2(e + 3e 2 ) 
dx 2 x dx 



1- 



2(3e 2 + 2c) 



= i • (13) 



At very large scales x *C 1 and at first order in the slow- 
roll parameters, the source term in Q setting S cm = 
S/k 2 reduces to 



8nGa 2 
~~k r ~ 

3 



x 2 p Lp 



nils' 

2(l + 2e-3e 2 ) 



2 [ H 2 



ft' -ft 



if 



n s 



3e dU s 



dx 



(14) 



All other contributions to the source term are suppressed 
by at least one factor x -C 1. The source is therefore 
completely dominated by the electromagnetic anisotropic 
stress Us- We evaluate lis m the next section; for this, 
we have to specify the generation mechanism for the elec- 
tromagnetic field which is operating during inflation. The 
above expression for the source includes also terms at first 
order in the slow-roll expansion: as will become clear in 
the following (c.f. Section HV| . for the problem at hand 
it is not enough to solve the Bardeen equation at lowest 
order in the slow-roll expansion, but we will need to go 
to first order. 

In the following, we will also need to solve for the cur- 
vature perturbation (see also (l9j). The curvature on 
comoving hypcrsurfaces £ is defined by 



C = * + 2(ft$ + #')/[3ft(l + w)} . 



(15) 



This variable has the advantage that it is known to be 
constant on superhorizon scales if the source is absent. 
Using the definition above, we can derive a first-order 
equation for £ which shows that even in the presence 
of the electromagnetic source, £ is conserved at lowest 
order in the large-scale expansion x <C 1, i.e., it does not 
contain a 1/x 2 term. Deriving (|15p . and eliminating $ 
via Eq. © and <£" with the help of the Bardeen Eq. ©, 
one obtains: 



c' = - 



2H 



3(1 + w) 



p v 



2ft + ^ 



Pern 
P<p 

3 \W 



foJ Pf a k 



,2 Hem j 



(16) 



where w = p v / ' p v = — 1 + 2e/3. The curvature perturba- 
tion is therefore sourced only at ncxt-to-leading order in 
the large-scale expansion i<1, i.e. it is of order x 2 ^/e. 
From the above equation, we see that the lowest-order 
solution for the Bardeen potential ^ is sufficient in order 
to calculate £ at next-to-leading order. 

It is also possible to derive an equation for £ which al- 
lows us to compute £ in a way independent of the Bardeen 
potential. This can be done directly from the Einstein 
equations in comoving gauge; the details of the deriva- 
tion are given in Appendix [21 and the resulting equation 
to lowest order in the slow roll parameters is 



dx 2 



2dC 
x dx 



+ c 



1 



ex 2 p v 



-6p e 



dp D 



dttf 



dx 



dx 



(17) 
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Comparing the source of the Bardeen Eq. (H~4l and the 
one of the above equation, we see that the former is 
by a factor x~ 2 larger than the latter. We therefore 
also expect the Bardeen potentials to be by a factor of 
about x~ 2 larger than the curvature perturbation, which 
is only sourced at next-to- leading order in a; <C 1. On 
the other hand, we note that the source of (|T7|) is larger 
in what concerns the slow-roll expansion: it is of order 
e _1 , while (fT4"|) is of order e°. One therefore needs to be 
very careful in dealing properly with the large-scale and 
slow- roll expansions, as will become clear in Section HVl 



III. THE SOURCE TERM OF THE BARDEEN 
EQUATION 

As discussed in the introduction, one of the simplest 
ways to generate an electromagnetic field by amplifica- 
tion of vacuum fluctuations during inflation is to break 
conformal invariance of the electromagnetic action by in- 
troducing a coupling between the electromagnetic field 
and the scalar field as 



S = - 



1 

167T 



. (18) 



with the Faraday tensor F^ = A Utfjl — A^ v , and A v the 
electromagnetic 4-vector potential. In the following, we 
adopt Coulomb gauge Ao(x.,rj) = 0, djAifcrj) = and 
follow the notation of Ref. [16[. From Maxwell's equa- 
tions, [■y/—gf 2 F flu ] tV — 0, we obtain an evolution equa- 
tion for the space components A4 (x, rf). In a cosmological 
background it reads [16| 



A'l + 2 ! -A\ - AAi 







(19) 



where A is the comoving spatial Laplacian. For a Fourier 
mode k, we simply have A = — fc 2 . The time evolution 
of the vector potential depends on the coupling function 
/(</?), and we adopt the following simple form for it @: 



m = h\ — 



(20) 



This choice is motivated on the one hand by simplicity, 
as it leads to simple power laws for the spectrum of the 
electromagnetic field. But it also includes the exponen- 
tial form first proposed by [BJ in the case of power law 
inflation, and it is sufficiently general to describe the case 
of generic single-field inflation in the slow-roll approxima- 
tion. We restrict to the values —2 < 7 < 2, which insures 
that the electromagnetic field remains subdominant and 
does not back react on the background expansion during 
inflation @, The value 7 = — 2 produces a scale- 

invariant (flat) spectrum for the magnetic field energy 
density, corresponding to a spectral index ns = — 3 for 
the magnetic field spectrum itself, as defined for example 
in Eq. (1.1) of [17| or Eq. flU below. 



With the above time evolution for the function /, the 
equation for the vector potential can be solved analyti- 
cally. Following [ill for the quantization of the electro- 
magnetic field, we expand the vector potential in terms 
of creation and annihilation operators b x (k) and b\(k) as 



d 3 k 
(2tt) 3 

b x (k)A(k, v )e 



E 

A=l 
ik-x 



el(k) 



(21) 



bl(k)A*(k,r,)e 



where ei(k), ea(k) arc unit vectors orthogonal to each 
other and to k, which represent the two polarizations 
of the electromagnetic field. It is convenient to define 
the new variable A = a(t]) f{rf) A(k,r}). Substituting 
Eq. ([201) and the expansion Eq. (J2J) into Eq. (fl9"]l. this 
latter can be solved in terms of the variable A as 



A{k,ri) 



Ci(7) J 7 _i/ 2 (x) + C 2 (7) J^ + i /2 {x) 



. (22) 

where x = \ki]\ = —krj, J„ denotes the Bessel function of 
order v, and C\, C2 are 7 dependent coefficients which 
are fixed as usual by imposing the initial condition that 
for subhorizon scales, — kr\ — > 00, the gauge field is in the 
Minkowski space vacuum [l6|. 

From the above solution for A we can infer the 
anisotropic stress IIs(k, rf), which appears in the source 
term in Eq. (JUJl. We have n s (k) = - 
3/2k l kjT c ^ ni , and the electromagnetic energy-momentum 
tensor is given by 



f 2 ( 

T v (an) — - I F F" 



1 = V El 



r- 



4ira 2 
(A ktj - A 3 . k ) 



AW 



(A k 



A,. 



(23) 



(V A A) 2 - A" 



Inserting the expansion Eq. ([2"Tj) into Eq. (f2"3"]l we find 



SV ' 11 2 fl 4 J (2n)\j^\ 3 



e Al (k')e v ,(k-k')/ 2 



Aik'^VfA^lk-k'l^V 



f 



f 



e A ^(k - k')(fcj - k' 3 ) - ev j(k - k')(k e - k' t ) 
e x t - e x i (kO^l A(k', V )A* ( |k - k' | , 77) 



6Uk')6t(k-k')+ c.c. 



(24) 



Here the 'c.c' stands for the three other terms with 
operators b x (k')b x ,(-k - k'), M k >v( k + k ') and 
b x (k')b\'(k' — k). More details are given in Appendix IB"1 
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From this expression for IT^k, 77) we can determine the 
source term in Eq. (1141) , which is a quantum operator act- 
ing on the electromagnetic vacuum. The source term is of 
the form S em (k, 77) = a x {k, r?)n s (k, r])+a 2 (k, rj)H' s (k, 77). 
How shall we proceed to compute the induced Bardeen 
potential? Naively one might simply want to use the 
vacuum expectation value of the operator S em as a clas- 
sical source term. However, this is not sufficient, since 
(0|ns|0) is independent of position and therefore does 
not contribute to the fluctuations. On the other hand, to 
work with the full fledged operator given in Eq. (|2~4")) is a 
bit unwieldy. The important point to remark, though, is 
that in order to solve the Bardeen Eq. © we only need 
to know the time dependence of Us (k, 77) , which allows 
us to calculate also (k, 77) and therefore the full source 
term. It turns out that, to determine the time depen- 
dence of ITs(k, 77), the easiest way is to first evaluate its 
power spectrum, 

(0|4(q, 7,)n s (k, 77)|0) = (2^) 3 P n (fc, r?M(q - k) , (25) 

where the <5(q — k) is a consequence of translation invari- 
ance and the spectrum P]i(k,r)) depends only on k = |k| 
due to the isotropy of the quantum vacuum. The details 



J 



of the calculation are given in Appendix [B] Here we only 
want to stress that, from Eq. (1241) . the quantum operator 
ns(k, 77) is of the form 



n s (k,77) =J2HKv)Oi(k) 



(26) 



where H^ are deterministic functions of time, and the 
operators Oi(k) do not depend on time (we formally 
perform the integral in d 3 k'). In the power spectrum 
(0|n^(q, 77)n,s(k, 77)|0), only one type of operators O, is 
such that (0|O](q)Oi(k)|0) 7^ 0, namely those which 
first generate two modes and then destroy them. The 
term which generates and destroys first a k' and then a 
q'-mode only contributes to the zero-mode, not to the 
fluctuation. There are two terms which give a nonzero 
contribution, and both give the same result. In Ap- 
pendix [B] it is shown that finally the anisotropic stress 
power spectrum can be written as the convolution of the 
magnetic, electric and Poynting vector power spectra (c.f. 



Eq. (TBIO 



P n (k lV ) = 18irf 4 f 
Jo 



> ^{a 1 ^)P E (k',r ] )P E (\k-k%r ] ) + a 2 (^P B (k',r 1 )P B (\k-k'\,r ] ) 



(2tt) s 



(27) 



+ cj 3 { 1 )P EB {k',r 1 )P EB {\k-k'\,r t ) 



Here we can neglect the contribution coming from k' > 
1/1 7/1, because for subhorizon modes the Bessel functions 
J v {k'\r]\) and J„(|(fc — fc')77|) which enter in the integrand 
(through Eqs. (|24p and (f2"2"|0 oscillate, and the result is 
damped. The prefactors 171(7) > a 2(l) and 173(7) depend 
somewhat on 7 but are always of order unity, and come 
from the angular integrals which cannot be evaluated an- 
alytically (see Appendix [Bjl . The power spectra are de- 
fined by 



results 



Pf 



Peb 



47T— -|^,77)| 2 , 
J a 

1 M(M) V 

a 4 \ f 

k ( A{k, n)\' 
fa 4 \ f 



4tt 



■A*(k, 77) 



(28) 
(29) 
(30) 



(0|fl i (q,7/)B;(M)|0) 
(0|£7,(q,77)£?(M)|0} 
(0|£ 4 (q,77)i?*(k,77)|0) 



M,)(2^) 3 P s (fc,77)<5(q-k) 
fc I fc,)(2^) 3 P B (fc,77),5(q-k) 



( S i3 

(Sij 

(\s l0 iki){2K? P EB {k,n)8{ci-K) 



Since we are interested in the solution for VE' of Eq. (TT3")) 
at large scales, supcrhorizon modes, we only need to com- 
pute the source for x = \kn\ < 1. We can then expand 
the Bessel functions in Eq. (|22|) for 1 C 1 and in this 
limit the solution becomes 



and have been calculated, e.g., in Refs. @, [l6| with the 



A(k,n) 



1 



ci(7)a; 7 + d 1 (j)x 



7+2 



Vk 

+ c 2 (7)2 a " 7 + d 2 (7^ 3 ~ 7 



(31) 



6 



with 



ci(t) 
did) 

C2(7) 

d 2 (7) 



a -i7T7/2 



cos(7r 7 ) 2(T-3)r(7 + l/2) 
ci(t) 



7+1/2' 

e i7T( 7 +l)/2 



cos(tt7) 2(7~T)r(3/2 - 7) ' 
22(7) 



3/2-7 



(32) 
(33) 
(34) 
(35) 



Depending on the value of 7, different terms dominate 
in the expansion (|3T|) . leading to different results for the 
magnetic, electric and Poynting vector spectra, and con- 
sequently also for Pjx(k,rj). For the power spectra we 
obtain, for x < 1 



P B {k,ri) 
Pe(M) 

PEB{k,1]) 



if 7 < 1/2 
|c 2 | 2 x 2 - 2 t if 7 > 1/2. 



T 2 ^ 

±E* J (WW^ 7+2 if 7 < -1A 

/ 2 a 4 }(l-2 7 ) 2 |c 2 | 2 a;-^ if 7 >-l/2 

47rfc 



(36) 



(37) 



-2|ci|-' a ,2'Y+l 



7 + 1/2 

l(27- l)c*c 2 

(27-l)|c 2 | 2 a; 1 - 2 ^ if 7 >l/2 



/ 2 a 



if 7 < -1/2 
if - l/2<7<l/2 



(38) 



As an example, we compute the spectrum of the 
anisotropic stress generated by the magnetic field Pb, 
i.e. the second term in the sum (|2T)) . For 7 < 1/2, we 
take the first line in Eq. (|3"6"|) . Using the general formula 
Eq. (|D2|) from Appendix [D] to approximate the convolu- 
tion, we obtain 



P 



n(B) 



9|ci| 4 a 2 
47r 2 |77| 5 a 8 

f 2 7 +l 

J (4+2 7 )(5+4 7 ) 

I 5+47 



(39) 



c 5+4 7 if _ 2 < 7 < -5/4, 
if - 5/4 <7< 1/2. 



For 7 > 1/2, we take the second line in Eq. (|36[) and we 
obtain 



n(B) 



9|c 2 | 4 g 2 1 
47r 2 |?7| 5 a 8 9-4 7 ' 



if 7 > 1/2 . 



(40) 



Similar expressions for the anisotropic stress generated 
by the electric field and the cross term are computed in 
Appendix ICl 



From these expressions, we can evaluate Pa(k,rf). 
Comparing the scaling of the magnetic, electric and 
cross term contributions, we can identify three different 
regimes. For —2 < 7 < —5/4, the magnetic field always 
dominates for x < 1, and we can neglect the electric 
contribution and the cross term. For —5/4 < 7 < 5/4 
all the contributions are of the same order of magni- 
tude. This follows from the fact that for these values 
of 7 the integrals over k' are dominated by the upper 
bound 1/1 7/|, leading to a white noise spectrum (see Ap- 
pendix [D]). Finally, for 5/4 < 7 < 2, the electric field 
contribution dominates and in principle we could ne- 
glect the magnetic field contribution and the cross term. 
However, at the end of inflation, when the Universe en- 
ters the radiation era, conductivity quickly becomes very 
high, meaning that the electric field decays rapidly (see 
for example |20(). The only remaining contribution to 
the anisotropic stress then is due to the magnetic field. 
Therefore, in this case we keep both the electric and the 
magnetic contribution in the anisotropic stress. Putting 
everything together we find the following power spectrum 
for the scalar anisotropic stress potential: 



J 



Pn(k,v) 



( |c 1 |V 2 (l+ 27 ) 5+ 4 7 
(4+2 7 )(5+4 7 ) x 



47r 2 |77| 5 a 8 



|4 

Cl| <T2 



IggllVl 



4|ci| (7 3 



if - 2 < 7 < -5/4 

if - 5/4 < 7 < -1/2 

if - 1/2 < 7 < 1/2 

if 1/2 < 7 < 5/4 



|ci| 4 <t 2 , |(l-2 7 )c 2 | 4 qi <t 3 |cic 2 | 2 (1- 27 ) 
5+4 7 ' 5-4 7 ~r 5 



\c 2 \ 4 <? 2 , |c 2 | 4 a-i(l- 27 ) 4 |c 2 | 4 g 3 (l- 27 ) 
9-4 7 ' 5-4 7 + 7-4 7 

|c 2 |V(l-27) 5 ^5-47 + |e 2 |V 2 



^ (4-2 7 ) (5-47) 



9-4 7 



if 5/4< 7 <2. 



(41) 



The expressions above diverge at the boundary of their validity. Except in the case \j\ = 2, this is simply be- 
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cause our approximation for the integrals derived in Ap- 
pendix [D] breaks down: the true integral would remain 
finite. For | — ^ | = 2 instead, this is the usual log divergence 
of a scale-invariant spectrum. Note that the magnetic 
field contribution for 7 > 5/4 (the term proportional 
to (T2 in the last line of the above equation) is smaller 
than the one from the electric field by a factor a; 47-5 , for 
x < 1. Hence the effect of this contribution to the dy- 
namical evolution of the Bardeen potential during infla- 
tion is negligible with respect to the effect of the electric 
field: we can therefore neglect it in the computation of 
the source during inflation. However, as previously men- 
tioned, only the magnetic field survives beyond inflation 
and generates the anisotropic stress, because of the high 
conductivity during the radiation era. We have therefore 
to take into account the magnetic field contribution for 
the solution of the Bardeen equation in the radiation era, 
see Section liVBI 

Furthermore, note that the prefactor |?7| _5 has the cor- 
rect dimension: since [Tls(x)] = and the Fourier 
transform is a volume integral, we have 



[t 



H^k)H s (k'))] = [<5(k-k')Pn] 



With [S(k - k')] = [l 3 ] this implies [P n ] = [£~ 5 ]. 

For the remainder of this paper, we use a simplified 
expression for the power spectrum Pn. For all values of 
7 e [—2,2], Eq. (|4"Tj) is of the following form, for super- 
horizon scales x < 1: 



cm. 



g.5-4171 if 5/4 < | 7 | < 2 



else, 



(42) 



where Cn(7) is a dimcnsionlcss parameter of order unity. 
This expression for Pn allows us to simplify considerably 
the analytical evaluation of the metric perturbations. It 
involves several approximations, but the spectral shape 
is correct and the evaluation of the prefactor is the best 
possible if one proceeds analytically. 

In order to compute the source of Eq. (fT4")) , we also 
need to determine Hg(k, 77); expression (|4"2")l helps us in 
this task. We note in fact that on supcrhorizon scales 
and at lowest order in the slow-roll expansion, the time 
dependence of the power spectrum of Hs(k, 77) is given 
by P n oc \-q\ 2a with 



4-2| 7 |, | 7 |>5/4 
3/2, | 7 |<5/4 

min{4-2|7|,3/2}. 



(43) 
(44) 



Moreover, from Eqs. (|2"4")l and (f3~Tj) one finds that on 
large scales x < 1, the operator n,g(k, 77), expressed in 
terms of the variables x and k, is simply a power law in 
x. Recalling Eq. ([2"6")l . we can simply write n,g(k, x) = 
xTa Ylt=i nj(k)Oj(k), where the entire time dependence 
is collected in the prefactor. Hence the power spectrum 
must go like x 2m : 

"'fiiWI^OlOftqJOiWlO), 

(45) 



(2^) 3 o(q-k)P n (fc) 



setting i = l for the only kind of operator which sur- 
vives. Eq. (|42p therefore implies m = a, so that we find 
the simple relation valid at lowest order in the slow-roll 
expansion: 



n's(M) 



■n s (k,77), atO(e°). (46) 



Note that at next order in slow-roll, n,s(k, 77) gets an 
extra time dependence as \ij\ 4e . We now have all the 
ingredients to determine the source term in Eq. (TT4"|) at 
lowest order in slow-roll: 



.7" 



3 , n s (k,77) 

-(2- a) 



(47) 



In Section llVDl to integrate the ^-equation (|T7| . we shall 
also need the spectrum of the electromagnetic energy 
density and its derivative. These can be obtained along 
the same lines as above, with the same qualitative results. 
We are not repeating the details. Defining as usual 

(OlpL(q^)Pcm(k, 77)|0) = (2^) 3 P om (fc,7?)<5(q-k), (48) 

and substituting Eqs. (f2"Tj) and (|2"2"j) in p om (x, 7/) = 
— T c ° m (x, 77), one obtains an expression anal- 
ogous to Eq. ([2"7]). where only the constants 
01(7), 02(7) and 173(7) are different. The final 
result can again be parametrized as 

Cem(7) ( x 5 - 4 ^ if5/4<| 7 |<2 
a 8|r?|5 \ 1 



P em (M)^^irK else ' (49) 



leading to 

Pem(k,?7) 



\V\ 



p em (k,77), atO(c°) (50) 



with a as defined in Eq. (|4"3"]l . 

With this, we can now proceed to solve Bardeen equa- 
tion. 



IV. RESOLUTION AND MATCHING 

A. Bardeen potentials during inflation 

We now want to compute "J at large scale by solving 
Eq. ([15]) with the source term given by Eq. (|47]l. The 
homogeneous solutions of Eq. ((T3)) arc "Ji = x v J v {x) and 
* 2 = x p J- v (x), with p = 1/2 + 0(e,e 2 ) and v = 1/2 + 
C(e, e 2 ) . The solution of the inhomogeneous equation can 
be computed using the Wronskian method, and reads 



*inh(k,x) 



2 sin(z/7r) 



dx 



X S'cm(k, X ) ^ ^ 



J v (x')J- v (x) - J v (x)J- v {x') 



(51) 



where x m = |fc77i n | is the initial time when the source 
starts to act, i.e. horizon exit xw. ~ 1, and it is such that 
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x ln 3> x. We denote the Bardeen potential by ^inh in 
order to indicate that this expression is a quantum oper- 
ator acting on the electromagnetic vacuum: we will then 
have to relate it to a stochastic variable after inflation 
in the usual way, see the discussion later in this section. 
Also in this case, the solution becomes highly squeezed 
after horizon exit. Only the dominant mode of the Bcsscl 
function, ^2 = x v J^ v (x) remains relevant. Since we are 
interested in the large-scale solution x -c x- m < 1, we 
can expand the Bessel functions in Eq. (|51|) for small ar- 
guments. Moreover, to integrate Eq. (|5ip we only need 
to know the time behavior of the source term, which is 
proportional to x a ~ A (c.f. Eq. (|47|l ). We obtain 



*inh(k,x) 



(2-a)(3-a) 
x 2 p v 



5cm (k, x) 



Note that, since < a < 3/2, the prefactor 



(3 -a) 



(52) 
(53) 

(54) 



is positive for all values of 7. Neglecting the decaying 
mode of the homogeneous solution ~$>i ~ s/^kx/A, we 
obtain the general large-scale solution for the Bardeen 
potentials during inflation, 



b(k) 
b(k) 



A Hi 

x 2 p v 
/3-3IIg 
x 2 p v 



and 



(55) 
(56) 



where b(k) is the usual inflationary solution at large 
scales, the homogeneous 'growing mode' (which is con- 
stant in time). To obtain Eq. (|56[) we have used Eq. ([7]) 
which yields 



JLIk 

x 2 p v 



(57) 



The above solutions are the sum of two uncorrelated 
quantum operators. The first, b, acting on the inflaton 
vacuum, and the second, proportional to n,g, acting on 
the electromagnetic vacuum. As mentioned before, these 
quantum variables must be identified with classical per- 
turbations having stochastic amplitudes by means of a 
quantum to classical transition: this is explained, for ex- 
ample, in [Ul^. For simplicity we identify these vari- 
ables with their r.m.s. amplitude, i.e. the square root 
of the volume factor in wave number space times their 
corresponding spectra: this is not a conventional choice 
but it allows to simplify considerably our formulas. From 
the corresponding power spectra, defined as 



(0|St(q,i 7 )S(k,i 7 )|0) 
|* t _(q,7 ? )#_(k,r,)|0) 



(2 7 r) 3 n(fc,7 ? ) ( 5(q-k), (58) 
(27r) 3 P*(fc, ?7 ) ( 5(q-k),(59) 



we therefore define the dimensionless metric perturba- 
tions: 



b^b= ^k 3 P b . 

An advantage of this somewhat unconventional definition 
is that both ^-(fc, 77) and b(k,rj) are dimensionless and 
provide a good measure for the fluctuation amplitude at 
comoving scale k. Analogously, for the anisotropic stress 
power spectrum we introduce the dimensionless ratio tt^j 
by 



[«5(M)] 

n~(k,ri) 



k 3 Pu (8ttG\ 
\'SH 2 J 

H 2 



k 3 P u 



Zm 2 p 



C n {l)x a . 



(60) 
(61) 



The superscript ~~ indicates that we evaluate the quantity 
in the inflationary era (as opposed to the radiation era, 
see Section llVBp . Note that here the Hubble parameter is 
taken at lowest order in slow-roll: H — H/a ~ 1/ (01/71), 
c.f. Eqs. (TIT)]) . With this definition, the power spectrum 
of the source term becomes 

(0|SL(q,r/)S om (k, „)| ) = (27r) 3 P s (A:, V )6(q - k) ,(62) 



(63) 



With Eqs. (|52l) and ((54)) . we find the relation among the 
power spectra 



k 3 R 



*i»h 



(2-a) 2 (3-a) : 



-k 3 P s 



(64) 



We can now rewrite the solutions Eqs. ([55]) and (|56| in 
terms of the classical variables, which we understand as 
their r.m.s. amplitudes, both for the inhomogeneous part 
of the solution and for the inflationary part of the solu- 
tion. We obtain with Eq. (l64l) 



(65) 
(66) 



#_(a;) ~ b(k)+/3^- and 



$_(x) ~ b(k) + (p-3) 



0; 



within the approximation that the magnetic field per- 
turbations and the inflaton perturbations are uncorre- 
lated, (0|* mh (fc)6(fc)|0) = (note that this is violated in 
second-order perturbation theory [H[). Therefore, the 
power spectra of and $_ are simply the sum of the 
inflationary power spectrum and the inhomogeneous one. 

The inhomogeneous part of the solutions Eqs. (|6"5)) and 
([66]). n^/x 2 , behaves like x a ^~ 2 . If 7 = -2, i.e. when 
the magnetic field energy density generated during in- 
flation has a scale-invariant spectrum, one has a = 
and therefore the inhomogeneous mode grows in time 
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like x~ 2 . One may wonder whether this leads to too 
large metric fluctuations, but it is not the case: consider- 
ing the ratio of the Weyl tensor C^ va p cx /c 2 (<£> + and 
the Ricci tensor, cx H 2 ~ rf 2 one finds for the ra- 
tio of typical components of the Weyl, respectively Ricci, 
tensor 1181 



Weyl 



Ricci 



inh 



inh J 



Cl n < 1 



(67) 



The last inequality is a consequence of the fact that we 
require the electromagnetic field to be subdominant dur- 
ing the inflationary era, so that II5 ~ p cm <C p v . 



B. Bardeen potentials during the radiation era 

After inflation and reheating, the Universe enters the 
radiation-dominated phase. Filled by a fully ionized 
plasma of relativistic particles, the Universe becomes 
conductive, in contrast to the inflationary phase during 
which there arc no free charges. The conductivity of the 
Universe is very high, so that the electric field disappears 
rapidly [f| [li| |20[ . Wavelengths of cosmological interest 
are much larger than the horizon scale at the end of in- 
flation. For them, the transition to the radiation era and 
the dissipation of the electric field can be considered as 
instantaneous. 

The evolution of the Bardeen potential in the radiation 
era has already been studied in detail in [TtJ and [Hj]. 
Einstein's equations can again be combined into a second- 
order equation for "P, that reads (see Eq. (B5) in (l7| ) 



(68) 



where again we identify the metric perturbations with 
their r.m.s. amplitudes, and fij^ = V k 3 Pu / p ra d is the di- 
mcnsionless magnetic anisotropic stress parameter in the 
radiation era. It is constant in time during the radiation 
era as both, the radiation density and B 2 scale as a~ 4 . 
Its value depends on 7 and is given by the magnetic field 
contribution of Eq. (|4"Tj) , i.e. by the part proportional to 
02 , evaluated at 77 = 77* . This magnetic part is indeed the 
only one that survives in the highly conductive radiation 
era. The general solution to Eq. (jSB"]) at large scales is 



3W 



and 



$+(x) = * 



*1 



(69) 



where ^0 an d ^1 are two arbitrary constants that need to 
be determined by matching the solutions in the radiation 
era to the one during inflation. 



C. Matching 

We have found the solutions for the metric potentials 
iff and $ in the presence of an electromagnetic field both 



during inflation (at lowest order in the slow-roll expan- 
sion) and in the radiation era. The solutions in the radi- 
ation era are known once the initial conditions are spec- 
ified. In order to have solutions valid through the whole 
evolution of the Universe, we need to match properly ^ 
and $ at the transition from inflation to the radiation 
era. 

The initial conditions in the radiation era are obtained 
by matching the solutions given in Eqs. (|65[) and (|66[) 
to Eqs. As usual, we match the classical potentials 

among themselves. We are only interested in wavelengths 
much larger that the duration of the transition. For these 
scales, the transition can be considered as instantaneous. 
Within this approximation, the equation of state param- 
eter w experiences a discontinuity at the transition, it 
goes from roughly w ~ —1 to w = 1/3. Furthermore, the 
electromagnetic field anisotropic stress is discontinuous 
when 7 > —5/4, due to the fact that the electric field 
contributes significantly during inflation for these values 
of 7, while it vanishes in the highly conductive plasma of 
the radiation era. 

In [T3| it has been shown that to match solutions 
through a discontinuity of the energy-momentum tensor, 
we have to impose that the induced 3-metric and the ex- 
trinsic curvature remain continuous on the spacelike hy- 
persurface of the transition S. We follow this procedure 
here. 

The most general metric containing only scalar pertur- 
bations is given by 



(70) 



ds 2 =a 2 | - (I + 2A)dif + 2B A drjdx 
+ [(I + 2C)5ij + 2E tij ]dx i dx j 



A convenient choice of coordinates to fix the match- 
ing conditions is to define the transition hypersurface 
through 77 = const. The time coordinate 77 is related to 
the original one by a gauge transformation 77 = T)+T [Til ]. 
On the {77 = const.} slices, the continuity of the induced 
3-metric and of the extrinsic curvature requires the con- 
tinuityof E, C, B — E' and & - HA through the transi- 
tion [14| . Considering the gauge transformation proper- 
ties of the different metric components [l8j], this implies 
in terms of the original perturbation variables 

[E}±=0, [C + HT} ± = 0, (71) 

[B — T — E']± = , [—HA + C- (H 2 - H')T]± = , 

where 



F ± = lim[F(+77* + e) - F(-rj* - e)} , 



(72) 



77* being the time of the transition: the end of inflation 
happens at 77 = —77* and the radiation phase is estab- 
lished at 77 = +77*. This also ensures the continuity of 
the Hubble parameter, the unperturbed extrinsic curva- 
ture: at lowest order in the slow-roll expansion one has 

Hini(r/ = -77*) = 7/7 1 = H rli d(ri = = K* ■ 
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Our original gauge is longitudinal gauge with E = B = 0, 
A = $ and C = With this Eqs. (O become 



[T]± = , [tf]± = , 

+ + (H 2 - H')T]± = 



(73) 



Inserting the solutions Eqs. (|S5|) and (|6"9")l into Eq. (|T3"|) . 

and using that {U 2 -U')- = tU\ and {U 2 -U') + = W.% 
one can determine the constants \&o an d in terms of 
T. Note that we set {H 2 - H')~ = Ml even though we 
are performing the matching at lowest order in slow-roll, 
since the order of T is as far not determined. A brief 
computation yields 



*o = b(k) + 



e - 2 



-n,r 



2/3 finO*) 



(74) 



— 



e-2 (3n^(x*) flt(x^) , N 

■^-H*T+^ nV 2 * y -3 m 2 ' . (75) 



If one inserts the above constants in the radiation so- 
lution Eqs. (|69[) . the following problem for the metric 
perturbations becomes manifest: since < a < 3/2, the 
term 



•ff* 



xT 2 > 



H !■ 



1/2 



which can be very large for small x», enters in the con- 
stant mode V&o of both and and contributes to 
the fluctuation amplitude as 

(*^* ) 



If this term is not 'compensated' in any way, it leads to 
very large fluctuations. Indeed, if one evaluates again the 
same ratio as in Eq. (|67[) but this time in the radiation 
era, one finds 



Weyl 



Ricci 



rad 



(76) 



which can become very large with the expansion of the 
Universe. However, we now show that this term is exactly 
compensated by the term proportional to T in Eq. ([74|. 
if we choose a proper hypcrsurphace for the transition 
and deal with the slow-roll expansion properly. 

To go on we need to specify a physically meaning- 
ful transition hypcrsurphace X, in order to determine 
T. Since inflation is driven by the scalar field, which 
controls the background evolution of the Universe, it is 
reasonable to assume that inflation ends when the en- 
ergy density of the scalar field reaches a certain value. 
In this case, the hypersurface of transition is the one 
of constant energy density of the scalar field, such that 
p (p (fj,^i) — p v + Sp,p = const, at constant fj. This requires 

to choose T = —Sp^/p'^, so that 5p v = [lij]- We can 
compute T explicitly at the end of inflation, for r/ = —77*. 



Using the scalar Einstein equations to eliminate the in- 
flaton perturbations we find 



n*T = 



-1 



3(1 



2$_ 



2*'_ 



2fc 2 *_ 



(77) 



Here = Pem/ P<p is the density parameter of the elec- 
tromagnetic field which is of the order of x 2 ^- like the 
last term in Eq. ([77)) : these two terms can therefore be 
neglected at lowest order in ir* <C 1. The first two terms 
in Eq. (|77[) are related to the curvature perturbation £ 
defined in Eq. (fT5"j) . Since \1/ and W do not jump, to low- 
est order in x* the jump in C, is related to the one in T 
by H*[T\± = - [C] ± . The condition [T}± = for the 
{pip = const.} hypersurface then implies 



0. 



(78) 



Hence the constant density matching requires that the 
curvature perturbation is continuous at lowest order in 
x* <C 1 (but it can be discontinuous at ncxt-to-leading 
order). Inserting Eqs. (|55|) and (|66|) into Eq. (|77|) . we find 
to lowest order in <C 1 and in the slow-roll expansion 



b(k) 



(79) 



It appears therefore that, if the hypersurface of transition 
is the one of constant energy density for the inflaton, T 
is of C(e _1 ) at lowest order in the slow- roll expansion. 
With this choice for T, the solution Eq. in the ra- 
diation era for <^ 1 is completely determined, with 
coefficients: 



*n = 



I(i-HW ) + ^_^ 

3 V ej 3 xi xi 



(80) 
(81) 



From these equations we see that the magnetic field 
anisotropic stress does not contribute to the leading order 
in slow-roll C(e _1 ) . The following order O(e ) cannot 
be trusted: this order would indeed receive contributions 
from T at next-to-leading order in slow-roll that we did 
not take into account in Eq. (|79[) . To be consistent it is 
therefore necessary to solve Bardeen equation up to next- 
to-leading order in the slow-roll expansion, as anticipated 
111 Section HU 

To solve Eq. (fi~3|) at next-to-leading order in the slow- 
roll expansion, we proceed as in Section IIV Al In the fol- 
lowing we present the solution and the matching for the 
inhomogeneous part only, sourced by the electromagnetic 
field; the relevant inflationary contribution in the radia- 
tion era is the one given in Eqs. (|80|) . ((8T|) at lowest order 
in slow-roll. For the inhomogeneous solution at next-to- 
leading order, we have to account for the fact that the 
source Eq. (|47p has now a time dependence of the form 
x a+2e . We find the inhomogeneous part of the solution 



(82) 
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Wc then identify the Bardccn potential quantum opera- 
tor with its r.m.s. amplitude as previously done, and we 
find the expressions (equivalent to Eqs. ([65]). (j66|) up to 
next-to- leading order 0(e)): 

*r nh (x) ~ (3 \l + (l + 2 log (i))cl § (83) 



£(84) 



*&h(*) * [P - 3 + (P + 2{fi - 3) log(£) 



The logarithmic factors arise because we have kept defi- 
nition (|61|l for £7^ , which is at lowest order in slow-roll. 

The matching proceeds as before with the difference 
that the continuity of the unperturbed extrinsic curva- 
ture % implies that the radiation phase is now established 
at time rj = 77* /(l + e): 



1 



■Hrad (77 



1 



H* . (85) 



Inserting solutions ([55"]) and (8"l]l into Eq. ((77)) . wc obtain 
for the electromagnetic contribution to the gauge trans- 
formation variable T at next-to-leading order: 



b(k) 



(86) 



Inserting this expression in Eq. (|74[) , we find that at next- 
to-lcading order, H*T cancels exactly the magnetic mode 

in * given by (2j8/3) Qn( x *)/ X l ( ci - Ec l- W)- so that 
in only the inflationary contribution remains. TTie 
matching conditions, Eqs. (|73[) . insure therefore that the 
magnetic mode oc J7^/a; 2 does nof transfer into the con- 
stant mode ^0 at the transition to the radiation era. Note 
that, in order to reach this conclusion, it is absolutely 
necessary to perform the matching. The decaying mode 
of the solution for ^ + in the radiation era which is pro- 
portional to f2jij does not give any information on vf^ (c.f . 
Eq. (M])). Therefore, the calculation done in Ref. [19[ 
and, in particular, their solution Eq. (33) are not suffi- 
cient to claim that no electromagnetic contamination to 
^0 is present after inflation. 

As a result of the matching, we see that the dangerous 
behavior inferred in ([76) is not present. Since the de- 
caying mode in Eqs. (|69[) can be neglected, we arrive at 
the result that the inflationary electromagnetic field con- 
tributes to the metric perturbations in the radiation era 
through ^0 only at next-to- leading order in 1 < 1, i.e. 
at order O(x ) instead of 0(x~ 2 ). Therefore, in order to 
find the relevant effect of the inflationary electromagnetic 
field, we would have to go to the next-to-leading order in 
the large-scale expansion. However, solving Bardeen Eq. 
(fTB"]) analytically at next-to- leading order in x = fc|?7| -C 1 
and in the slow-roll expansion is highly nontrivial. The 
electromagnetic contribution to the metric perturbations 
can be found more easily by solving for the curvature 
perturbation £ which, as shown at the end of Section [TTl 
is only sourced at next-to- leading order ini< 1. 



D. The curvature perturbation 

The metric perturbations in the radiation era are given 
in Eqs. (69)) : neglecting the decaying mode, the relevant 
contribution generated by the electromagnetic field is in 
the constant mode ^0, but only at order O(x ) as demon- 
strated in the previous section. The simplest way to ob- 
tain ^0 is to solve for the curvature perturbation f . In- 
deed, inserting Eqs. (fBT)]) into the definition of ( Eq. ([15)) 
one finds that 



3^o 
2 



(87) 



Since £ is continuous at the transition from inflation to 
the radiation era, as demonstrated in Eq. (|78[) , this means 
that ^0 is given by the value of the curvature at the end 
of inflation 



3 



1 + 7 ' 



where wc denote by C* the contribution sourced by the 
electromagnetic field during inflation, which we need to 
determine. 

Note that the solutions Eqs. (|6"5"|) and (jBB"]) do not allow 
to compute £*• Indeed, inserting them into Eq. (|15|) we 
find that the term HQ + ty' vanishes. Naively one might 
think that this implies simply £_ = and that ( has a 
contribution at order 0(x~ 2 ) ■ O(e ); however, this would 
be in contradiction with Eq. (fTB")) , which shows that £ is 
only sourced at next-to- leading order in x <C 1. In reality, 
this contradiction is solved if one inserts in definition (|15[) 
the Bardeen potentials at order e and not simply e , be- 
cause the term "H$ + ty' enters with a prefactor of order 
0(l/(w + 1)) = Cie' 1 ) in (15]). Inserting solutions (83]) 
and (|84|) one verifies that the order 0(x~ 2 ) ■ O(e ) van- 
ishes, and we expect that the same thing happens at any 
following order in e (note that Eq. (fTB")) is indeed valid at 
any order in slow-roll). Consequently, the relevant con- 
tribution of the electromagnetic field is of the order 
O(x ) ■ 0(e~ 1 ), as anticipated at the end of Section ITT1 
To evaluate (* using Eq. (fT5"]) . we would therefore need 
to go to the next order in and <&_: we evade this 
by computing directly. This can be achieved either by 
solving Eq. (]16[) and inserting solution ([65]) for or by 
solving directly Eq. (TJJ . 

Let us start by solving Eq. (fT7|) . The source term is 



St = 



1 



ex 2 p v 



-6p e 



dp e 



dllf 



dx 



dx 



(89) 



involving not only the anisotropic stress but also the elec- 
tromagnetic energy density and its derivative, which are 
given in Eqs. (|49l50p . The electromagnetic energy den- 
sity and the anisotropic stress have the same dependence 
on time and wavenumber, i.e. they go as x a . The in- 
tegration of Eq. ([TTj) becomes then straightforward: the 
details of the solution are given in Appendix [A] Identi- 
fying again the quantum operators C and p cm with their 
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r.m.s. amplitude as done in Section IIV A[ one finally 
obtains the solution 



H 2 1 



9mp e 



{a - 6) 2 C 2 m + 2a(a - 6)C p n + a 2 Cg 

— log (ar/xi n ) if a = 
xf n /a if a ^ 



1/2 



(90) 



where X[ n ~ 1 denotes horizon exit (when the source 
starts to act). We have introduced the coefficient C p u 
denoting the amplitude of the cross term arising from 
the correlation (0|/?J m IIs|0): the cross-correlation has the 
same spectral dependence, the electromagnetic 

energy density and anisotropic stress. We can define, 
analogously to Eq. (|6ip. 



^em(M) 



l >2 r 



H 2 
3m 2 p 



C em ( 7 ) x a 



(91) 



It is interesting to note that in the scale-invariant case, 
a = 0, the anisotropic stress does not contribute to £r h 
to lowest order in the slow-roll expansion. Furthermore, 
only in the scale-invariant case there is a logarithmic 
build up of £in h , whereas in all other cases C;„h * s con " 
stant on large scales: it is generated at horizon crossing 
and then stops growing. The Bardeen potentials ^_ and 
on the other hand, keep growing outside the horizon 
and soon become large: this renders longitudinal gauge 
badly adapted to the problem at hand. Conversely, in the 
comoving gauge, for example, all perturbation variables 
remain small during inflation (c.f. Appendix |A")) . 

We could have solved for fr, by means of Eq. (JTS]) and 
inserting solution (|55|) for In this case, the source 
term in (|16|) contains also a contribution form the Poynt- 
ing vector. However, this can be eliminated in favor of 
Pom and II5 using the momentum conservation Eq. (IA10[) 
derived in Appendix [XJ and the two approaches give fi- 
nally the same result for Ci„h- 

Neglecting the decaying mode and using so- 

lution with C* = CinhC 2 -*)' tne Bardeen poten- 

tials in the radiation era are completely determined (c.f. 
Eq. ra): 



2 / 1 

2 / 1 

3 1 + 7 



b(k) 
b(k) 



2C* 
3 

3 



*S and 



(92) 



Comparing these solutions with the ones obtained at 
large scales from a magnetic field generated causally dur- 
ing a primordial phase transition [l7j . we see that the 
only difference is the constant mode 



2C* _ ^yi 

3 \ mp 



log (x*) if a = 
1/q if a ^ 



(93) 



(neglecting factors of order one). This constant mode is 
generated by the curvature perturbation at the end of 
inflation. In the case of a primordial phase transition, 
the curvature perturbation has no contribution from the 
electromagnetic field before the phase transition. The 
continuity of the curvature at the phase transition forces 
therefore = 0. During inflation, however, the cur- 
vature is dynamically generated by the electromagnetic 
source and at the end of inflation it transfers into the 
constant mode V^Oj adding a new term to the Bardeen 
potentials. 

In solutions (|92p . the mode sourced by the magnetic 
field anisotropic stress in the radiation era fij^ would 
dominate at large scales. However, once neutrinos start 
free-streaming, this term is exactly compensated by the 
neutrino anisotropic stress, and leaves no visible effect in 
the CMB, see 17 . The observationally relevant contribu- 
tion from iljij is therefore given by the so-called 'passive' 
mode, arising at next-to- leading order in the large-scale 
expansion [231 ] - However, the passive mode is present 
also for a magnetic field generated by a causal process 
such as a phase transition; in the case of a magnetic 
field generated during inflation, one has to add the term 
proportional to £* to the passive mode. The full solu- 
tion for the Bardeen potential in the matter era from an 
inflationary magnetic field becomes therefore (in this dis- 
cussion we neglect the contribution of the compensated 
mode, which is anyway subdominant with respect to the 
passive mode): 

» +0 . >*j-|(i +7>« + x-f SJ i ** (S). 

(94) 

where x cq = kr] cq with ?y eq denoting the time of equality, 
and x v = kr\ v with r\ v denoting the time of neutrino de- 
coupling. The last term is the passive mode taken from 
Eq. (6.10) of [It}, where we have set for the time of cre- 
ation of the magnetic field 77 — 77^ . Let us compare the 
amplitudes of the new constant mode 2(*/3 and of the 
usual passive one. For sufficiently red spectra, a < 3/2, 
we can safely set ~ f2jj and equally Q~ m ~ Q+ m 
(c.f. discussion at the end of Section riVB[) . The ampli- 
tude of the constant mode is therefore proportional to 
(H*/mp) 2 ~ Q^/x", while the one of the passive mode 
is proportional to fi^ ~ f2~ m . Therefore, for a scale- 
invariant electromagnetic field with a = 0, the difference 
among the amplitude of the two modes is only due to the 
presence of the logarithms and of 1/e (c.f. Eq. ([93])). The 
logarithm appearing in the term 2£*/3 corresponds to the 
number of c-folds between horizon exit of a given scale 
k ~ 1 /r]i n and the end of inflation, while the one appear- 
ing in the passive mode corresponds to the time interval 
between the end of inflation and neutrino free-streaming. 
If, on the other hand, ^ a < 3/2, the amplitude of the 
passive mode is significantly suppressed with respect to 
the constant inflationary one due to the extra factor x" . 

Note also that the new term 2£*/3 is of the same 
order of magnitude as the inflationary power spectrum 
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Pb — (i?*/mp) 2 /e, up to the large log in the scale in- 
variant case a = 0. However, coming from the square of 
a Gaussian field, it is genuinely non-Gaussian and may 
lead to a significant /nl, see Refs. [T3. l24j|. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have computed the impact of pri- 
mordial magnetic fields generated during inflation on the 
Bardeen potentials in the subsequent radiation era. We 
have solved the Bardeen equation both during inflation 
and in the radiation era. We have used the inflation- 
ary solution as initial condition for the radiation era, 
by matching the solutions at the transition in such a 
way that the induced metric and the extrinsic curva- 
ture on the p = constant hypcrsurface remain contin- 
uous through the transition. This matching procedure 
uniquely determines the solution after inflation. 

We have shown that the leading-order contribution to 
the Bardeen potential at large scales, \E'_ ~ Q,^/ (k-q^) 2 , 
is transferred entirely to the decaying mode in the radia- 
tion era. Consequently, the inflationary electromagnetic 
field contributes only at next-to-leading order, i.e. at or- 
der O(x ). A similar situation can appear in bouncing 
universes, like e.g. the ekpyrotic universe [25|, [26|: dur- 
ing the contracting phase, the growing mode can become 
very large. Only if this large mode is entirely transferred 
to the decaying mode of the expanding Universe, do per- 
turbations remain small; this full transfer can be achieved 
by the continuity of the curvature perturbation £. Here 
the situation is different, because the large mode is due 
to a source. But we have found that even in this case, £ 
remains small and its continuity is sufficient to guarantee 
that the large, superhorizon mode does not contribute to 
the constant ^o- This shows that a gauge which exhibits 
large metric fluctuations, as the longitudinal one, can be 
very misleading: in the transition to the radiation era 
the Bardeen potential <f> jumps by a huge amount, while 
£ remains continuous. 

Performing the calculation at order O(x ), we have 
found that the Bardeen potential in the radiation era 
contains a constant term proportional to the amplitude of 
the electromagnetic energy density at the end of inflation, 
which takes the form 



(fc?7*) Q e 



mp 



(95) 



If the electromagnetic field generated during inflation is 
scale-invariant, this contribution is enhanced logarithmi- 
cally by a factor of about — log(fcr7*). Note that the ori- 
gin of this constant term in is quite different from 
the case of a magnetic field generated by a causal pro- 
cess, like, for example, at a phase transition in the early 
Universe. In this latter case, in fact, the matching condi- 
tions insure that the electromagnetic field perturbations 
are compensated on superhorizon scales. 



In addition to the above contribution to there is 
the one from the anisotropic stress of the magnetic field 
during the radiation era oc Sl^ / x 2 , which is however com- 
pensated later on by neutrino free-streaming. As shown 
in [I?) (see also [I?]]), the magnetic field anisotropic stress 
is not compensated as long as the neutrinos are still cou- 
pled to radiation and do not free-stream. However, when 
neutrinos decouple, they develop an anisotropic stress 
that counter-balances the one of the magnetic field, re- 
moving the magnetic mode at large scales 28j. Only 



the metric perturbations at next-to-leading order remain, 
which correspond to the so-called 'passive' mode [23| . 
This compensation also takes place if the magnetic field 
is generated during inflation. However, as we have shown 
with this work, in addition to the mode which later is 
compensated, oc £l^j/x 2 , and to the passive mode, infla- 
tionary magnetic fields also lead to a constant mode as 
in Eq. ([9"5]l. 

Let us estimate the amplitude of the constant mode 
for the most interesting case of a scale-invariant mag- 
netic field. In this case, the anisotropic stress does not 
contribute to the result since a = 0. The electromagnetic 
energy density is 



(96) 



so that T* 



Furthermore, H 2 ~ T^/m 2 p < 

rap (f^) 1 / 4 . The largest observationally allowed mag- 
netic fields have an amplitude B ~ 10~ 9 G, leading to 
n- m ~ l(T 7 and T* - 10" mp. For the conformal 
time at the transition we find ry* = H.~ x = (a*7J*) _1 = 
T*VOem/(2bmp). Setting fco = Hq ~ /(rapy/TE^), 
we obtain for the present Hubble scale ko = Hq 
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-28 . 



so that — log(x»o) ~ 64. Therefore, the amplitude of 
the metric perturbation, vp™ 1 ~ — Q~ m log(x*)/e, at the 
wavenumber corresponding to the present Hubble scale 
is enhanced by the number of e-folds of inflation after the 
present Hubble scale has exited the horizon. 

To estimate the effect of an inflationary magnetic field 
on the CMB at large scales, we can set very roughly 
AT/T ~ the full observationally relevant contribu- 
tion is then given by the sum of the passive and the con- 
stant modes, as given in Eq. (|94[) . In the scale-invariant 
case a = 0, the logarithmic enhancement of the passive 
mode is —3/5 log(T*/T„) ~ —27, and the one of the new 
inflationary contribution is of the order of the number of 
c- folds, therefore in total larger by about a factor 2/e. 
We obtain then that the amplitude of the effect on the 
CMB is increased with respect to the naively expected 
amplitude ~ Q~ m by nearly 2 orders of magnitude due to 
the large logarithms. In other words, it might be possible 
to detect inflationary magnetic fields in the CMB down 
to about 10 _10 G instead of the usual limit of 10 _9 G. 
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The limit of about 10 _9 G was obtained, for example, in 
Refs. [1^] which, in order to be model-independent, only 
evaluate the effect of the compensated mode, correspond- 
ing to "P™ 1 ~ &em- Furthermore, contrary to the inflaton 
perturbations, the contribution from the magnetic field 
is inherently non- Gaussian as it is quadratic in the field 
amplitude. This leads to a nontrivial bispcctrum, sec 

Refs. dsim. 

It remains to be investigated whether it is consistent 
to ignore the gauge symmetry breaking which is neces- 
sary for the proposed mechanism to work, or whether it 
leads to other consequences which have to be taken into 
account. For example the effects of a possible longitudi- 
nal gauge mode have to be studied. However, apart from 
fine-tuned examples like the one proposed in Ref. [3(| , at 
present this inflationary scenario seems to be the only vi- 
able candidate for primordial large-scale magnetic fields. 
If a causal generation mechanism is to be successful, a 
strong inverse cascade is needed in order to move corre- 
lations from small to larger scales. This inverse cascade 
has to be more efficient than the one proposed e.g. in (3lj . 
which has been shown to be insufficient [32j . The issue 
of the evolution of helical magnetic fields is still unsolved 
and interesting research in this direction is ongoing, see 
e.g. 0. 



In this gauge the perturbed Einstein equations read 



HA-H' L -*f = 



m 2 A - SHH'r. - k 
AmGa 2 



(A2) 
(A3) 



ikP_ 



k 2 A + H'+ + 2HH' T -k 2 \H L 



em j Hcmj 



8ttGci 2 I1 s 



HA' + (2H> + U 2 )A -^A-^.f HL + ^ 



(A4) 

ZHH'r 



-H'{ = ^Ga 2 l- 



<p> 2 A 



Pcm ik J / , \ 

+ ~ + ~akl 1 \ 3nqcm J + q ^Jj 

(A5) 



These equations can be rewritten in terms of the curva- 
ture perturbation £ which is given by 



(A6) 
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Appendix A: An evolution equation for the 
curvature perturbation 

1. In comoving gauge 

Here we derive an equation for the curvature pertur- 
bation £ in the comoving gauge. Since the curvature is 
gauge invariant, we can use the solution to compute the 
constant mode of the Bardccn potential as shown 
in the main text. We work in the gauge comoving with 
the total fluid, i.e. such that q\ ot = qf m + qf = and 
Bi = 0. The metric is (in Fourier space, we use the 
notation of fl8^ 



ds 2 =a 2 \ - (l + 2A)dr) 2 



[(1 + 2H L )8ij + 2H T Yij] dx l dx 3 



(Al) 



Using Eq. (|A2[) to express A in terms of £ and Eq. (|A6|) 
to write Ht in term of Hl and £, one finds the following 
set of coupled equations for £ and Hl'- 



- 2U 



H 



C - mH' L + k 2 c 



(A7) 



AirGa 2 { -p, 



ok 1 



cm j Hem j 



3C" 



m - — C' - k 2 C + ZH'L + 6nH' L = -8TrGa 2 U s 



H 



C~^C + H'l + 2UH' L = 



(A8) 
(A9) 



emj + Qemj 



Combining Eq. (| A8|l with Eq. (|A9[) one finds a conserva- 
tion equation for the magnetic field 



Pern 



2n.c 



(A10) 



Deriving Eq. (|A7|) and combining it with Eqs. (|A8[) 
and (|A10[) gives a second-order evolution equation for 
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the curvature £ 

w 




(All) 



At lowest order in slow-roll this equation becomes 

^xow'xj.V 8nGa 2 f p' em U s ' 

C + 2UQ + k C = 7^— 1 6p om + -y- + -y- 

(A12) 

The source term in Eq. (|A12[) has no contribution propor- 
tional to Hs/k 2 . Hence as expected it remains small at 
all scales providing that the magnetic field is small. The 
curvature seems therefore a more suitable variable than 
the Bardeen potentials to describe perturbations in the 
presence of a primordial magnetic field during inflation. 
In terms of the variable x = — krj Eq. (|A12[) reads 



d 2 ( _ 2 d( + £ 
dx 2 x dx 



1 



ex 2 p v 



dp G 



dx 



St 



dx 
(A13) 



As we argue in the main text, the time and k dependence 
of p em are the same as for Us, i.e. oc x a . We can therefore 
rewrite p eirL and II5 simply as 



Pern 
P v 



A em x a and 



Ik 

p v 



(A14) 



where A om and An denote the amplitude of the quantum 
operators which do not contain any x dependence and 
are therefore irrelevant for the integration of Eq. (|A13I) ■ 
With this, the source term for the curvature becomes 



St 



„a-2 



(a - 6)A em + aA n 
The homogeneous solution to Eq. (|A13I) is 

C = cyx 3/2 J 3/a (x) + c 2 x 3/2 J_ 3/2 (a;) . 



(A15) 



(A16) 



The solution of the inhomogeneous equation can be com- 
puted using the Wronskian method. Integrating from 
x- m ~ 1 to < x <C 1 one finds 



Cinh 



1 



(a-6)A cm + aA n 



log (x/x- m ) if a = 
xfja if a ^ ' 
(A17) 



From the above solution, it appears that when one identi- 
fies the quantum operator £; n h with its r.m.s. amplitude, 
in general there will be also a cross term due to the cor- 
relation among the electromagnetic energy density and 
the anisotropic stress, (0|Aj m A n |0) (c.f. Eq. of the 
main text). 

2. In synchronous gauge 

The ^-equation of motion can also be obtained in syn- 
chronous gauge. There the metric is 

ds 2 = a 2 {-dr/ 2 + [(1 + 2H L )5 ij + 2H T Y ij ]dx i dx j } . 

(A18) 

And the perturbed Einstein equations read 



H' L + ^ - 4ttG ( 



iak 3 



~ Qemj 



H r 



H, 



4nG [ip' Sip' + V v a 2 Sip + a 2 p em ] 



2 



2UHL — k 2 [ Ht, H — ) = 87rGa 2 n 



H 



4vrG 



ip'^Sip' — V. v a 2 Sip + a p e 



2 , 
- 3 a 2 U s 



(A19) 
(A20) 

(A21) 
(A22) 



Combining these equations to eliminate Sip and using 
that in this gauge 



<~jb(* + x)-(* + x> 



one can again derive Eq. (|A13[) for £ in the slow-roll ap- 
proximation. 



Appendix B: Computation of the anisotropic stress 

Here we present a detailed computation of the 
anisotropic stress power spectrum. First we need to com- 
pute the anisotropic stress ns(k, rf). As an example, we 
focus on one specific contribution, namely 



IIs(x, 77) 



1 P 



2 87m 4 



A' i (x,r ] )A , i (x,r ) ) + 



(Bl) 



Using the expansion in Eq. ([2~Tj) for Ai , the Fourier trans- 
form of the anisotropic stress reads 
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n s (M) 



2a 4 



d 3 xe lkx 



d 3 k 



d 3 q 



(27r) 3 j (2?r) 3 J^e^OOe^Cq) 



0(^)'e^ + 6l(k')(^)' ( 



-ik -x 



Mq) 



M(g,q)V 



e iq ' x + &l(q; 



-jqx 



(B2) 



This equation contains four different products of creation and annihilation operators. For each product the integral 
over d 3 x can be performed, e.g. 



&A(k>v(q) / d 3 xe iX ( k+k ' +q) = (2^) 3 ( 5(k + k' + q)6 A (k')6A'(-k-k'), 
and similarly for the three other terms. With this the anisotropic stress becomes 

___ f 2 f d 3 k' /^(fc',7 ? )V 

17 v ; AA' = 1 



(B3) 



V / 



'e A „(-k-V)^<' k + k 'l'')V 



& A (k')6l,(k + k')+e Vi (k'-k) 



& A (k')& A ,(k-k') 



/ 



M(ik'-kU)V 



6 A (k')&A'(-k-kO 



f 



6|(k')&v(k'-k) 



(B4) 



Performing the same calculation for the other contributions in Ilg we obtain Eq. (|2~i|) . 

We can then compute the spectrum of the anisotropic stress (0|llt(q, 77)II,s(k, ry)|0). The only operators that 
contribute to the spectrum arc bb^bb^ and bbb^b^ . The first operator corresponds to the case where one mode is 
created and destroyed and then a second mode is created an destroyed, and the second operator corresponds to the 
case where two modes are created and then destroyed. The contribution from the first operator reads 



(0|(6 Q (q')6^(q + q')) t &A(k')^(k + k')|0) ^ S xx ,6 aa ,S(k)5(q) . 



(B5) 



This term does only contribute to the zero-mode k = q = and it has therefore no effect on the fluctuations, i.e. on 
the Bardeen potential. The contribution from the second operator is 

(0|& a ,(q-q')&a(q')&i(k')^(k-k')|0) 

= (2ir) 3 5 aX 5(q' - k')(0|Mq - q')&A<(k - k')|0) + (0|M<1 - q')& A (k> a (q')& A ,(k - k')|0) (B6) 
= (2nfS aX 5 a ,yS(q! - k')S(k - k' + q' - q) + (2ir) 6 5 aX ,6 a , x 5(k' + q' - q)5(q' + k' k) , 

where for the first equality we have used the commutation relation between b a and b\. The two terms in Eq. (|B6[) 
give both the same contribution to the power spectrum, which becomes 



(0|n s T (q,ry)n s (k, ?7 )|0) = ^(q-k) / d 3 k' { B( (k, k')/ 
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fA(\k-k'\,r,)\' 
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Bi{k,k')\A(k',r))\ ■ U(|k-k'|,77)| +2B 1 (k,k')B 2 (k,k')A(k',ri)A(\k-k'\,ri) 



x/ 2 



A*(k', V )\' /-A"(|k - k'|,7 7 )\ ' 



/ 



where 
^i(k,k' 



E 

A,A' = 1 
A,A' = 1 



f 



tik? - ^ ) 



k'k 3 



e xl {k')k[-e Xl {k')k[ ■ e A ^(k-k')(* J --*5)-evi(k-k / )(fc<-&J) 



(B7) 



(B8) 
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Equation (|B7[) contains integrals over the direction of k' that are difficult to compute exactly. However, since we 
are mainly interested in the scaling of the anisotropic stress with rj and k, but not in its precise numerical value, we 
approximate these integrals with 



1 

-in 



dn k /£?(k, k') 



#-kj" )V 



0-1(7) 



(B9) 



with (71(7) a constant which depends somewhat on the index 7 through the solution for A but which is always of 
order one since the terms in B\ contain only unit vectors. Similarly we approximate 



i- J dn k , J B 2 2 (k,k , )|^(|k-k'U)| 2 ~a 2 ( 7 )fe' 2 |fe-fe'| 2 \A{\k-k'\,r,)\ 



(BIO) 



and 



A J d^B^k^B^k^Adk-klr,) ( ' 4 * (|k j klv) J ^ a 3 (j)k'\k - k'\A(\k - k'\, V )\ f , 

(Bll) 

with (T2 (7) and 03 (7) two constants of order unity. With these approximations only the integral over the wavenumber 
k! remains and we find for the anisotropic stress power spectrum 



9 -47T f^k'^dk'j 



(2tt) 3 



■A(fc',r?)V 
/ 



A(\k-k'\,r,)\' 



f 



+ a 2 { 1 )k' 2 \k - k'\ 2 \A{k' , v )\ 



x \A(\k - k'\,r])\ + cr 3 (7)fc'|fc - k'\A(k',r])A(\k - fc'U)/ 2 



^(fc',r,)V /A*(\k-k'\, V )\' 



(B12) 



/ 



where we neglect the contributions coming from k' > l/\r}\ because for these the Bessel functions in A, J v (k'\r)\) and 
J v {\{k — k')r)\), lead to oscillations and the result is damped. This equation can be written in terms of the power 
spectra for B, E and EB, that have been calculated, for example, in Refs. (5, 
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(B14) 
(B15) 



With this we find 



Pn(k,v) -18^/ 4 / 
Jo 



/•Vkl k'9iik' i 

j { a x {i)P E {k! t ri)P E {\k - fc'U) + a 2 (j)P B (k', v )P B (\k - k'\,rj) 



(2tt)s 



(B16) 



+ <T3h)PEB(k',Tl)PEB(\k-k , \,Tl) 



The calculation of the electromagnetic energy density 
spectrum P cm {k,r]) is analogous. The only difference is 
in the angular dependence of the individual terms which 
yields different parameters ci, a 2 and (73. 



Appendix C: The various contributions to the 
anisotropic stress power spectrum 



Here we compute the spectrum of the anisotropic stress 
generated by the electric field Pe and the cross term Peb ■ 
Let us start by the electric part. For 7 < —1/2, we use 
the first line in Eq. Q37p . Approximating the convolution 
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with Eq. (|D2|) in Appendix [D] we find 

9|2 Cl |V 1 



P, 



n(-E) 



4tt 2 (7 + l/2) 4 |?7| 5 a 8 9 + 4 7 



(CI) 



For 7 > —1/2 we take the second line in Eq. (|57)) which 
gives 



Pn(* 



9|c 2 |V(l-2 7 ) 4 
47r 2 |ry| 5 a 8 



5-4 7 

if - 1/2 < 7 < 5/4 



1-27 



(4-2 7 )(5-4 7 ) 

_ if 5/4 < 7 < 2 



„5-4 7 



(C2) 

The computation of the cross term involves three cases. 
For 7 < —1/2, we use the first line in Eq. ([38]) which 
gives 



Pi 



9MV3 



27+2 



(5+2 7 )(7+4 7 ) x 

if - 2 < 7 < -7/4 



47r 2 (7 + l/2) 2 |?7| 5 a 8 

7+4 7 

if - 7/4 < 7 < -1/2. 

(C3) 

For —1/2 < 7 < 1/2, we take the second line in Eq. 
and wc obtain 



P 



9|ci| 2 |c 2 | 2 cT3(l-27)- 
47r 2 5|ry| 5 a 8 



(C4) 



Finally for 1/2 < 7 < 2 we use the third line in Eq. (|58"]) 
which gives 



Appendix D: Convolution integrals 

In the convolution f^ /M dk'k' 2 P x {k')P Y (\k - k'\) we 
usually have to integrate power laws. Hence these inte- 
grals are of the form 



I(a,p)(k) = 



i/M 



dk'k' a \k~k'f 



(Dl) 



Integrals of this type are very common when dealing with 
primordial magnetic fields; the standard way of approx- 
imating them has been given first in [Hj|. More refined 
analytical evaluations of integrals of this type are beyond 
the scope of this paper, see e.g. (28[. Wc first note that 
these integrals require a + 1 > in order to avoid an 
infrared singularity at k' — > 0. We then split the integral 
in its part < k' < k and k < k' < l/\rj\. We use that 
x = k\r]\ < 1, hence < k < In the first interval 

we approximate |fe — fc'| ~ k while in the second interval 
we set \k — k'\ ~ k'. With this approximation, which 
is certainly crude but retains the main characteristics of 
the behavior, we obtain 



7(a,/3)(fc)~ 
+1 a+8+lV n J 



a + 



1 



if a + P + 1 > 



_ I Q+/3 + 1 



(D2) 



9|e 2 | 4 a 3 (l-27) s 
4 7 r 2 |? ? | 5 a 8 



7-4 7 

if 1/2 < 7 < 7/4 



2-27 



(5-27) (7-4 7 ) 
^ if 7/4 < 7 < 2 



,7-4 7 



(C5) 



The result for Pjj(b) is given in Eq. ((40 



For the last ~ we have set x = k\rj\ and we use x < 1 
to determine the dominant contribution. Interestingly, 
such a convolution always either has a red spectrum, oc 
fe™, n = a + P + 1 < or it is white noise, oc fc°. Blue 
spectra cannot be generated by a convolution. If small 
scales dominate, the integral is dominated by the upper 
cutoff which yields a white noise behavior. 
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